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Abstract 

The effect of weak anisotropic (momentum-dependent) impurity scattering in 
unconventional superconductors has been investigated. It is shown that the 
anisotropic scattering can lead either to a small reduction or a small enhance- 
ment of the isotropic pair-breaking effect. The influence of the anisotropy 
of the scattering potential becomes significant for the order parameters with 
large Fermi surface average values. In that case an unexpected enhancement 
(up to ~ 10%) of the critical temperature Tc over the critical temperature in 
the absence of impurities Teg is predicted for a small impurity concentration. 
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There has been a great interest in studying the properties of unconventional supercon- 
ductivity in last decade. These superconductors are characterized by order parameter of 
lower symmetry than the conventional one, which only breaks U(l) gauge symmetry. Such 
order parameters are involved in almost all theories of heavy fermion and high temperature 
superconductivity.EF0 As the order parameter A (k) is a nontrivial real function of k its 
Fermi surface average (A (k)) strongly depends on the shape of the Fermi sheet and may 
take values between -1 and 1 under the normalization condition (A^ (k)) = 1. Most of the 
theoretical effort has been focussed on the case of (A (k)) = however, and only recently the 
unconventional superconductivity with a nonvanishing Fermi surface average of the order 
parameter has been discussed.Hl^ In the presence of impurities particularly, these supercon- 
ductors were suggested to show some interesting phenomena in the density of states, which 
should detectable in the thermodynamic and transport properties.i^lHllI 

In this paper, we go beyond the approximation of the isotropic impurity scattering and 
consider the weak anisotropic (momentum-dependent) impurity potential, which can be 
treated as a small perturbartion of the isotropic scattering. We show that a really interesting 
case of weak anisotropic scattering is when the Fermi surface average of the order parameter 
is nonzero. Particularly for large values of (A (k)) we obtain a surprizing increase of the 
critical temperature Tc with the impurity scattering rate in the range of small impurity 
concentration. 

Recent experimentalEUll as well as theoreticafl studies suggest the importance of the 
anisotropic impurity scattering issue in the interpretation of the Tc suppression due to im- 
purities in YBCO. As the anisotropy of the impurity potential is preasumably strong in 
this compund@'0 the use of a perturbation method may not be applicable. We believe, 
however, that the results of the present study may apply to some other anisotropic super- 
conductors like heavy fermions, for instance, where a highly complex Fermi surface^l (FS) 
and anisotropic order parameter may lead to a nonzero (A (k)).lii0 We take h = ks = 1 
throughout the paper. 
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We assume randomly distributed nonmagnetic impurities in an anisotropic supercon- 
ductor. Treating the electron-impurity scattering within first-order Born approximation 
and neglecting the impurity-impurity interaction,lii the normal and anomalous temperature 
Green's functions averaged over the impurity positions are given by 

^(..,1.) = ^^±i^ (1) 



A(k) 

C^2+^,2^|A(k)|2 



^(^.k)= ^_^ , (2) 



where the renormalized Matsubara frequency uj (k) and the renormalized order parameter 
A (k) read 

(k) = a; + im \w{k- k') {u, k') 3 (3) 

J (27r) 

r d^k' 
A (k) = A (k) +nj\w{k- k!) k') (4) 

J (277) 

In above u) = 7rT(2n + 1) (T is the temperature and n is an integer number), is the 
quasiparticle energy, rii is impurity (defect) concentration, w {k — k') is momentum depen- 
dent impurity potential and A (k) is the orbital part of a sing let§ superconducting order 
parameter defined as 

A(k) = Ae(k) (5) 

where e (k) is a real basis function of an one dimensional (ID) irreducible representation of 
a crystal point group or a linear combination of the basis functions of ID representations. 
We normalize e(k) by taking (e^) = 1, where < ... >= Jpg dSkn{k) (...) denotes the average 
over the Fermi surface (FS), n (k) is the angle resolved FS density of states normalized to 
unity, i.e. JpgdSkn (k) = 1, and Jp^^dSk stands for the integration over the Fermi surface. 



The impurity scattering potential is assumed to be separable and given by 



4 



\w{k-k')\^ = \w\^h(k)h{k') (6) 

with 

hik) = l + gik), |(7(k)|«l (7) 

which means that g (k) function represents a small anisotropic correction to the isotropic 
scattering potential. The above functions are normalized by taking (h) = 1 {{g) = 0). 
We note from Eq. and the form of impurity potential (Eqs. and (|^)) that uj 

is k-dependent. This means that the electron self-energy due to impurity scattering and 
consequently the quasiparticle life-time are anisotropic and change over the Fermi surface. 
Further, it yields from Eqs. and (0) that the impurity scattering may change the 

symmetry of the renormalized order parameter A (k) depending on the g (k) symmetry. In 
that respect our approximation differs from that by Markowitz and Kadanoff as well as 
that by Millis et al.0 where the authors assumed only a change of a degree of order param- 
eter anisotropy but not the anisotropy function itself. 



To proceed further, we restrict the wave vectors of the electron self-energy and pairing 
potential to the Fermi surface and replace / d^k/ {2tt)'^ by Nq JpgdSkn (k) / d^k, where Nq 
is the overall density of states at the Fermi surface. Using Eqs. (|I]), (^, (^, and in Eqs. 
(0) and (^) and performing the integration over (particle-hole symmetry of quasiparticle 
spectrum is assumed) we write 



(D(k) = uj[l + u{iu)h{k)] 
A (k) = A [e (k) + e{u)h (k)] 

where u (a;) and e (uj) functions are determined by the self-consistent equations 

l + u{Lu)h (k) 



u{uj) = T f dSkU (k) h (k) 
Jfs 



CD2+ |A(k) |2 



1/2 



e{u) = T I dSkU (k) h (k) — 

JFS 



e (k) + e {uo) h (k) 



cD2+ |A(k) 



1/2 



(8) 
(9) 

(10) 
(11) 
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and r = TrNoTiilwl'^ is the impurity scattering rate. The gap function is given by the self- 
consistent equation 

A(k') 



A(k) = -Tvy y(k,k')— , ~ 



(12) 



where V (k, k') is the phenomenological pair potential taken as 



y(k,k') = -roe(k)e(k') (13) 
Following standard procedure,ii we obtain the equation for the critical temperature Tc as 



<x)>0 



(14) 



with 



(/M)..o-tc^^.n(k)^ 



A(k) ' 
A 



(15) 



A=0 



where T^^ is the critical temperature in the absence of impurities and c^o (k) = uj (k)^^Q. 
Using Eqs. (^-(^) and neglecting terms of order (k) and higher, we obtain for A 



ijJq (k) = cij + r (1 + (yf (k)) sign (uS) 



(16) 



A(k) 
A 



e(k)+ (e) + (e^) 1- 



r 



J A=0 



\uj\ + r 



r 



\UJ\ 



[1 + 9 (k)] 



(17) 



We note from Eq. (|T^) that the impurity scattering induces anisotropy in the renormalized 
order parameter A (k) other than that of A (k) which is described by a function e (k) (Eq. 
(^)). There are two different sources of this additional anisotropy in A (k). First of them is 
a nonzero value of (e) which is determined by the symmetry of the order parameter and that 
of the Fermi surface, and the second reflected in Eq. (|1^) by g (k) function and coupling 
coefficient l^eg) is the anisotropy of the impurity scattering potential. 



Based on Eqs. (P^D, (|TB|), and (0) we obtain from Eq. (|T^) 
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(18) 



where tjj (z) and ^jj^^^ (z) are digamma and trigamma functions respectively. The first term 
on the righthand side of Eq. (|TH]) represents the isotropic impurity scattering effect on 
TciUni and the influence of the weak anisotropic scattering is reflected by the second term. 
Compared to the isotropic pair-breaking, the weak anisotropic impurity potential can lead 
either to a small (of order g) additional decrease of Tc when the term (2 (e) (eg) — {e^g)) is 
negative or to a small reduction of Tc suppression for a positive value of (2 (e) {eg) — {e^g)). 

We consider the case of a positive value of this parameter in more detail. For the sake of 
convenience we assume {e^g) = and then deal with the order parameter-impurity potential 
coupling coefficient (eg), which does not reduce the generality of approach but simplifies the 
notation. One would obtain the same results without the above assumption using the term 
(2 (e) {eg) — {e^g)) as a parameter with the values assigned to 2 (e) {eg) coefficient in this 
paper. The actual value of {e^g) depends on two factors: the symmetry of both e (k), g (k) 
functions and the shape of the Fermi surface. Thus the estimation of {e^g) parameter is 
difficult in the real systems, however, its value may vanish in some cases (see Appendix). 

We present the critical temperature T^. normalized by the critical temperature in the 
absence of impurities Tcq as a function of normalized impurity scattering rate T /2t{Tcq for 
different values of the impurity coupling coefficient {eg) in Figs. 1-3. As we are considering 
a weak anisotropic potential and have neglected terms of order {eg)"^ while writing Eq. (|T8|), 
we take values of {eg) upto ~ 0.3. We notice from these figures that the effect of anisotropic 
scattering is very weak for small Fermi surface average values of the order parameter (e) and 
may not be distinguished from the isotropic impurity scattering ( Figs. 1-2), nevertheless 
the influence of the anisotropy in the scattering potential increases with the increase in (e). 
For the large values of the Fermi surface averaged order parameter that is (e) ~ 0.9 (Fig. 
3) and for the coupling coefficient {eg) larger than approximately 0.1 we observe a very 
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interesting feature of an initial enhancement of Tc for small impurity concentration. This 
initial enhancement of over T^.^ may be understood by writing Eq. (p!8|) for small impurity 
concentration (i.e. r/27rTcQ ^1). We have in that case 



positive a and decreases with the impurity concentration for a negative a. This criterion 
may serve for estimation of {eg) value which gives rise to the enhancement of the critical 
temperature with scattering rate. For example for (e) = 0.8 {{e^g) = 0) this coefficient 
is (eg) ~ 0.23, and for (e) = 0.9 {{e'^g) = 0) the coefficient (eg) = 0.2 leads to the 
enhancement of about 10% for r/27rTc„ ~ 0.12 (as in Fig. 3). The required large (e) value 
suggests a substantial s-wave component in the order parameter. 

Although in above we have taken (e'^g) = 0, we expect, that even when (e^g) coefficient 
has a positive value, the factor a can be made positive by appropriately large value of (e). 
It is worth mentioning that the above results cannot be obtained with the assumption of the 
proportionality of the impurity potential anisotropy function g (k) to the order parameter 
orbital function e (k), since the constraint {g) = would yield (e) = 0. 

In summary, we have investigated the effect of weak anisotropic impurity scattering in 
unconventional superconductors. It is shown that the anisotropy of the impurity potential 
can lead to either a small increase or a small decrease in suppression due to isotropic 
scattering. The influence of weak anisotropic scattering becomes signiflcant for the order 
parameters with large Fermi surface average values (e). In this case the critical temperature 
is increased up to 10 per cent over T^q for small impurity concentration. This anisotropic 
scattering-induced Tc enhancement in unconventional superconductors with large (e) values 
is a novel feature of the present study. 

This work was supported by the Natural Sciences and Engineering Research Council of 




(19) 



where a 



(e)^ — 1 + 2 (e) (eg) — (e'^g). Therefore Tc is an increasing function of F for a 
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Canada. 



APPENDIX: 

We assume e (k) and g (k) as the basis functions of ID irreducible representations of the 
crystal point group. Therefore for a symmetry operation from this group Si 

5,^(k)=ft^(k) (Al) 

where gi is a number. Further, if there is a part P of FS which reproduces the whole Fermi 
sheet under the crystal group symmetry operations (FS=X^j SiP) and the average value of 
g (k) over this segment of FS does not vanish, i.e. {g)p 7^ 0, then from the normalization 
(g) — and from the relation 

i i 

we get 

Y.9i = (A3) 

i 

Taking into account that (k) is an identity representation for a real e(k), which means 
that Sic'^ (k) — e} (k), we calculate {e^g) as follows 

^-'9)-i:{S.(f9)),-^e'9)X^. (A4) 

i i 

Then Eqs. (A3) and (A4) yield {e^g) = 0. 

On the other hand, if the order parameter is given by a linear combination of a s-wave 
component and a basis function / (k) of a nonidentity ID irreducible representation of the 
crystal point group, that is 

e (k) = (e) + / (k) (A5) 

then {eg) = (fg) and under the same assumptions about g (k) we get (e^g) = 2 (e) (eg) 
which leads to a cancelation of the anisotropy of the scattering potential in Tc equation 
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(Eq. (|T8D). We conclude from this, that in order to observe the effect of weak anisotropic 
scattering, e (k) and g (k) must belong to Aig representation with g (k) normalized to fulfill 

{g) = 0. 

The above considerations, however, are based on the assumption that g (k) transforms 
according to the symmetry of the crystal point group and are not valid in the case of g (k) 
described by a symmetry other than that of the crystal lattice. 
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FIG. 1. Normalized critical temperature Tc/Tc„ as a function of the normalized scattering rate 
V/2'kTcq for (e) = 0.1. Prom the bottom to the top, curves are plotted for [eg) =0 (isotropic 
scattering), 0.2, and 0.3; {e^g) = is assumed. 
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FIG. 2. Normalized critical temperature Tc/Tc^ as a function of the normalized scattering rate 
r/27rrcQ for (e) = 0.5. From the bottom to the top, curves are plotted for {eg) =0 (isotropic 
scattering), 0.1, 0.2, and 0.3; (e^g) = is assumed. 
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FIG. 3. Normalized critical temperature Tc/Tcq as a function of the normalized scattering rate 
r/27r for (e) = 0.9. From the bottom to the top, curves are plotted for {eg) =0 (isotropic 
scattering), 0.1, 0.15, and 0.2; {e^g) = is assumed. 
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